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Abstract. Two problems about the faintest Small Magellanic Cloud (SMC) Cepheids are addressed. On one hand 
evolutionary tracks fail to cross the Cepheid Instability Strip for the highest magnitudes (i.e. I-mag~ 17) where 
Cepheids are observed; Mass-Luminosity relations (ML) obtained from evolutionary tracks disagree with Mass- 
Luminosity relations derived from observations. We find that the above failures concern models built with standard 
input physics as well as with non-standard ones. The present work suggests that towards highest magnitudes, 
Cepheids stars undergo a selection effect caused by evolution: only the most metal poor stars cross the Instability 
Strip during the "blue loop" phase and are therefore the only ones which can be observed at low luminosity. This 
solution enables us to reproduce the shape of the lower part of the Instability Strip and improves the agreement 
between observed and theoretical ML-relations. Some issues are discussed, among them Beat Cepheids results 
argue strongly in favor of our hypothesis. 
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1. Introduction 

Cepheids are variable stars located in the Color- 
Magnitude diagram (CM-diagram ) within the Instability 
Strip (IS) where pulsation phenomena take place via the 
mechanism. Cepheid masses approximately range be- 
tween ~ 3 M and ~ 15 M©. During the past decade, 
microlensing experiments as MACHO 1 , EROS 2 , MOA 3 
or OGLE 4 have produced a huge flow of data. As by- 
products of these observational programs, a large num- 
ber of new variable stars, and among them, Cepheids 
have been detected. OGLE 2 data provide a large and 
high quality sample of Cepheids belonging to the Small 
Magellanic Cloud (SMC) which we consider in this work. 
This extended and homogeneous data set has already put 
several shortcomings of the SMC Cepheid modeling into 
light: 

(1) the evolutionary tracks built with standard input 
physics and for a chemical composition Zo = 0.004 (Yq = 
0.251) as usually assumed for the SMC fail to reproduce 
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the observed Cepheid position within the CM-diagram for 
highest magnitudes (i.e. ~ 17 mag); 

(2) the Mass-Luminosity relation (M L puls hereafter) 
deriv ed from pulsation properties (see iBeaulieu et alJ 
l200ll) and ML-relation from evolutionary tracks (ML evo1 
hereafter) do not agree. 

The first problem arises because theoretical "blue 
loops" do not cross the observed IS over the whole Cepheid 
mass range. Indeed an evolutionary track for a mass about 
5 Mq crosses observed IS three times, the first time ("first 
crossing") is the faster one -e.g. ~ 0.01 Myr for an 5 
M Q model-, the second time ("second crossing") is the 
slower -e.g. ~ 0.20 Myr- and the third time ("third cross- 
ing") remains short, about ~ 0.01 Myr. These time scale 
considerations tell us that the majority of the observed ob- 
jects should be in the second crossing stage. During this 
phase, the star burns He in its inner regions. Third and 
second crossing both belong to the so-called "blue loop" 
excursion towards the blue side of Hertzsprung Russell 
(HR) diagram. Consequently theoretical blue loops should 
cross the observed IS for the entire Cepheid mass range. 
As we confirm in the first part of this work, theoretical 
tracks with standard physics and free parameters varied 
in a reasonable range are not able to provide blue loops 
which reach the observed SMC Cepheids at low magni- 
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tude, i.e. they are not able to model the low luminosity 
Ccpheids in the SMC case. 

The s econd problem has b een underlined -among 
others- bv lBeaulieu et alJ l|200l|) who found a strong dis- 
agreement between the Mass-Luminosity relation ML puls 
and ML ™ 1 in the SMC case. They determine a ML puls for 
the LMC and SMC using pulsation calculations, indepen- 
dently from evolutionary calculations: for a given Cepheid 
the mass M* and the luminosity L* are found iteratively 
solving an equation of the type P i thoo (M*, L+) = P° bs , 
where P° bs is the observed period (i — for fundamental 
pulsators and i = 1 first overtone ones) and P* heo the the- 
or etical one calculated w ith a pulsation code. The results 
of iBeaulieu et all ((200]Jl are based on calculations that 
assume a metallicity content Zq = 0.004 which is usually 
assumed to represent the mean metallicity of the SMC. 

We are therefore led, in the second part of this paper, 
to propose another possibility and show that it can rec- 
oncile both issues: the blue loops at low mass and the 
Mass-Luminosity relation problem. In our hypothesis no 
high magnitude (i.e. with a mass around ~ 3 M Q ) SMC 
Cepheids with a metallicity as high as the mean SMC 
value can exist because the evolution does not bring these 
stars far enough on the blue side to cross the instability 
strip. The observed high magnitude SMC Cepheids must 
therefore be undermetallic (i.e. Z ~ 0.001) with respect 
to the mean SMC metallicity (i.e. Z ~ 0.004). 

In Sect. |3 we recall the physical inputs used in our 
standard models, which are similar to what is found in the 
recent literature. We next compare our evolutionary tracks 
with OGLE 2 observed Cepheids within a CM-diagram. 
Follow ing: a method similar to that used bv lBeaulieu et alJ 
l)200l|) . we also compare ML ™ 1 and ML puls . In both cases 
we confirm the discrepancy. 

In Sect. |2| we discuss the above issues in view of the 
uncertainties of the standard models and discuss the effect 
of non-standard physics in cases when models including 
such physics are available. In Sect. 0] we compare models 
calculated with Zq — 0.001 with observations. Sect. [S] is 
devoted to discussion about the possibility of the existence 
of SMC low luminosity Cepheids with metallicity as low 
as Z Q = 0.001. 

2. Standard models versus observations 

2.1. Inputs for Standard evolutionary models 

Our evolutionary models are built with the ID Heny ey 
type code CESAM 5 originally written bv lMorell (|l997l) in 

which we brought several improve ments. 

The equation of state is from lEggleton et al ] lll973l) 
and the external boundary condition is defined in a 
simplified model atmosphere involving the Eddington 
T(t) law. The nuclear n etwork involve s 30 nu clear reac- 
tions, we have followed ISchaller et all ill 99^ who used 
the same network as in iMaederl ()l983|) for H-burning 

5 CESAM : Code d'Evolution Stellaire Adaptatif et 
Modulaire 



and iMaeder fc Mevnetl l|l987j) He-burning network sup- 
plemented with the 17 O(q,n) 20 Ne reacti o n. Nu clear reac- 
tion rates are from ICaughlan fc Fowler! lll988j) excepted 
12 C(a, 7 ) 16 0, 17 Ofo.^ ^F from ICaughlan et a,lJ |j 985h 
and 17 0(p,a) 14 N from lLandre et alJ l|l990h. More re- 
cent n uclear rates do exist: NACRE bv lAngulo et alJ 
l|l999|) . however adopted rate for 12 C(a,7) 16 is quite 
simila r to NACRE one ( a facto r of about two higher 
than ICaughlan fc Fowler] (|l988l) and about 80% of 
ICaughlan fc Fowlerl (|l988() one. 

The adopted mean chemical composition for the SMC 
is taken as X = 0.745, Y = 0.251 and Z = 0.004, corre- 
sponding to a met al to helium enrichme nt of AYq / AZo = 
2 (see for instance iPeimbert et al.l (12000^ ~) . to a primordial 
helium Y P = 0.243 Jlzotov et al1ll997l) and to [Fe/H] = 
-0.68 JLucket alJ Il998h. Elemental abundances corre- 



spond to the lcrevesse fc Noels! l)l993j) (GN93) mixture 
consistent with OPAL96 calculations. 

Ilglesias fc Rogers! 



Opacities are from Igl esias fc Rogers! (199E 
(O PAL96) for high temperatur es (T > 10,000 K) 
and I Alexander fc Ferguson! l|l994[) for cooler domains. We 
stress that the central chemical composition during the 
He burning phase differs strongly from GN93 (e.g. 50 % 
of 12 C and 50 % of 16 0). Thus we have use opacity tables 
allowing a variable composition in 12 C and 16 O with 
the aim of modeling the core as rea listically as p o ssible . 
These tables have been built with iMagee et~a ill 995h 
elemental opacities (Los Alamos). 

The convective flux is computed accor ding to the pre- 
scrip tion of the Mixing Length Theory l|Bohm-Vitensel 
1958). The mixing length value Zmlt -derived from solar 
calibration- is equal to 1.6 H p . We used Schwarzschild's 
criterion to decide if the energy transport is radiative or 
convective, and an extra mixing zone is added above the 
convective core (i.e. overshooting). The extension of this 
zone is taken to be l over = 0.2 H p (l over = a over H p ). 

For the transformation of theoretical quantities, 
(Mboz, T c ff) into absolute magnitudes and colors, we 
used the Bas el Stella r Library (PaSeL, version 2.2) of 
iLeieune et alJ 1^998) which provides color-calibrated 
theoretical flux distributions for the largest possible range 
of fundamental stellar parameters, T c g (2000 K to 50,000 
K), log.g (-1.0 to 5.5), and [Fe/H] (-5.0 to +1.0). 



2.2. Comparison between models and observations: the 
"Blue Loop problem" 



Cepheids data are from lUdalski et al ] lll999h . We have 
chosen to work with (V — I) colors for which more data 
are available. Fundamental and first overtone Cepheids are 
plotted in Fig. [TfbJ. A mean (V — I) reddening is taken 
from lUdalskfjl998h . E(V - I) = 0.08. The SMC distance 



modulus is fixed at 18.94 from lLanev fc Stobiel l|l994|) with 
an i nternal error of 0.04 mag, this is a well accepted value, 
e.g. lCroenewegenl l|2000h found 19.11±0.11 or 19.04±0.17 
depending on the photometric band. 
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Fig. 1. (a) Theoretical HR diagram showing our standard evolutionary tracks with masses 3.0, 4.0, 5.0 and 6.0 M Q . 
The adopted chemical composition is Xq = 0.745 Yq = 0.251 Zq — 0.004. The heavy elements mixture composing Zo 
is from GN'93. The core overshooting extension is 0.2 H p along the entire tracks, (b) CM-diagram showing OGLE2 
data (black dots) including the sample of Cepheids (fundamental and first overtone pulsators). Only 10 % of the 
stars belonging to th e SC 1 OGLE 2 field main sequence are displayed for (V-I)< 0.2. The adopted reddening is 
E(V - I) = 0.08 from llJdalski et all l)l999h . The cross on the right side indicates an estimation of errors: 30% of 
E(V — I) = 0.08 on color and 0.1 mag on I, for Cepheids. The effective temperature which should be reached by 
models with masses of about 3 is indicated by a vertical dashed straight line. The open diamonds are the sample 
of Cepheid s used bylLuck et al.l l)l998(l in their study of chemical composition, values for I-magnitude and (V — I) are 
those from lLuck et, al.l l|l998|) 's Table. 3 In panel (b), the evolutionary tracks of (a) are also plotted (solid lines). The 
straight dashed line A represents the limit below which the evolutionary tracks fail to model the observed Cepheids. 



In both figures (Fig. [Ja) and Fig. ^b)) evolution- 
ary tracks involving standard input physics are dis- 
played. FigQJa) is a HR-diagram showing log L/Lq versus 
logT c ff . A segment of a line shows the temperature which 
should be reached -according to OGLE observations- by 
the evolutionary tracks for a stellar mass of about 3 M ; 
we will mention this mark in further discussions. An es- 
timation of the uncertainties is also plotted in Fig. ^b): 
we estimate the error on the colors to be roughly 30% 
of E(V — I) (which is the typical variation of reddening 
within the sample) and assess an error of 0.1 dex on I- 
magnitude which roughly represents the distance modulus 
uncertainty. 

In Fig. [lib) we have also plotted the stars from 
iLuck et "all 1^998) for information. 

The general characteristics of these theoretical dia- 
grams are simil ar to those shown by se veral groups like 
Geneva one, see lCharbonnel et al.l 1 1993h . This is not sur- 
prising because these authors have used similar physical 
inputs. For instance, the effective temperature at the tip 
of our 3 M Q blue loop is logT B ff = 3 .728 which com- 
pares well with lCharbonnel et all jl99$\ logT c ff = 3.734. 
In all cases, the logT ff value is far from the required one 
of about 3.82, i.e. a temperature hotter by ~ 1200 K. 



The main features put in evidence in Fig ^ arc: 

(1) the main sequence position seems to be reasonably 
well reproduced by the models 

(2) the position of the blue tip of the 3 M blue loop 
is too red. For an I-magnitude corresponding to a mass 
of about 3 M Q , we can clearly see a bulge of Cepheids. 
In fact 93% of fundamental pulsators and 81% of first 
overtone pulsators are located between I m j n ~ 16.5 and 
I mol ~ 17.7. Such a large amount of objects -statistically 
significant- cannot be explained solely by ~ 4.0 M Q first 
crossing models. Indeed, for a 4.0 M standard model, 
the time spent during the first crossing is Tfi rst = 4.6 x 
10~ 2 Myr while the time it takes for the second and third 
crossing is T" seconc } cr ossing+ r third crossing = 3433 M y r - 
Hence, blue loops should cross the entire observational IS 
for the lowest masses. The adopted value of distance mod- 
ulus /i does not affect this conclusion. Indeed, e ven if we 
take e xtreme evaluations: [i\ = 18.66 ±0.16 fromlUdalskil 
l|l998|) and fi 2 = 19.05 ± 0.13 from iKovacsl l|2000|) . the 
evolutionary track for ~ 3 M Q does not extend through 
the observational IS. 
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Fig. 2. Mass-Luminosity relation derived from OGLE ob- 
servations for Z a = 0.004 corresponding to [Fe/H] ~ —0.7. 
Circles: fundamental pulsators, squares: first overtone pul- 
sators. Solid line: mass-luminosity relation from our evolu- 
tionary code with Zq = 0.004 and an overshooting amount 
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0.2 H n ; da shed line: mass-luminosity relation from 



iBono et al.l l|200i 

2.3. Comparison between observed and calculated 
Mass-Luminosity relations 



2.3.1. Deriving Mass-Luminosity Relations from 
observations 

In order to derive a ML-relation from the observations, we 
use a method very similar to the one used by Beaulieu et 
al. (2001). For each object we solve iteratively the equa- 
tion: 



Pt hco (M+,L+,T cB ,Y ,Z ) = P' l 



obs 



(i) 



where P° hs is the observed period value (i — for 
fundamental pulsators and i — 1 for first overtones) 
and P i theo (M*,i*,r eff ,y ,^o) the theoretical one, com- 
puted with the Florida LNA 6 pulsat ion code which is 
a Castor type code (see ICastorlll97l|) . During the iter- 
ative process, M* is adjusted in order to match P° hs 



and p* hco ; for a given iteration M* is fixed and we solve 
the following set of equations where the unknowns are 
(log T off , log i?*,log L*): 



logT eff = 3.9224 + 0.0046 log g 
- 0.2470 (V - I - {R v - 



- 0.0012 [Fe/H] (2) 
Rj)E{B - V)) 



2.5 log U = Msmc - V + Rv E{B - V ) + BC + 4.75 (3) 



4waRlT* B 



r M * 



(4) 



(5) 



Eq. Incomes from lKovacsl l|2000Fl Eq.2, in which we brought 
absorpt i on cor rections. Eq.[3]is the Beaulie u et al.'s Eq.2. 
iKovacsl l|200d|) who made interpolations of ICastelli et alJ 
stellar atmospheres models to convert magnitudes 
into bolometric and effective temperature into colors. The 
luminosity L+ is in solar units, 



BC = 0.0411 + 2.0727AT - 0.0274 log 5 



0.0482[Fe/H] 
-8.0634AT 2 



Linear Non-Adiabatic 



and AT = logT c ff — 3.772. The magnitude and color V, 
(V— I) are from OGLE observati ons, msmo has been take n 
equal to 18.9 consisten t ly wi th lLanev fc Stobid l)l994|) . 
Following ludalski et alJ l|l999|) we took E(B - V) ~ 0.08, 
R v = 3.24 and i?/ = 1.96. 

In order to apply this method one has to select the 
data. Indeed on the CCD detector, a Cepheid may be 
"blended" with another star, the magnitude of the object 
being shifted towards lower magnitudes. These "over- 
luminous" objects lead to wrong couples (log M*, log T*), 
therefore it is crucial to reject from the sample the stars 
suspected to be blended with other object(s). From 
OGLE data we have extracted amplitudes of pulsation in 
B, V and I bands and then derived Amplitude — logP re- 
lations. The criterion to suspect that an object is blended 
is the following: if a given object has a magnitude lower 
than the mean magnitude (at least 0.2 mag lower) given 
by the Magnitude — logP law and an amplitude lower 
than the mean amplitude given by Amplitude — log P 
relation, this object is rejected. Moreover we have also 
rejected some objects which appear to be suspiciously to 
red. Finally we retain 1177 fundamental pulsators and 
709 first ove rtone pulsators and ob tain similar samples of 
objects than iBeaulieu et al 1 (I2001J) . 

Fig- El displays the resulting ML-relation derived with 
an assumed metallicity Zq = 0.004. W e did not find 
signifi cative differences with results from IBeaulieu et alJ 
1120011) . 

2.3.2. Uncertainties in Derived ML-relations 

For years the question of Magellanic Clouds distance 
has been a subject of debate. There were supporters 
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for "short" distance scales -e.g. IStanek et all l|l998|) with 



ulmc = 18.065 ± 0.031 ± 0.09 mag 


- and for "long" dis- 


tance scales -e.g. 


jancv & Stobie 


1994 with m.iun = 


18.53 ± 0.04 mag- 


Cioni et alJ 


(2000) derived a distance 



mod ulus for the LMC Ut.mc? = 18.53 ± 0.04 ± 0.08 mag 
and iMould et al.1 l|200o|) (HST Key Project Team) have 
adopted /ilmc = 18.50 ± 0.04 ± 0.15 mag. We made 
a test with "short" distance (i.e. /ismc = 18-7 mag); 
this leads to differences in mass of 5 log M/Mq ~ —0.1 
and in luminosity 8logL/ Lq ~ —0.08 (consistently with 
Eq.|3|. We dismiss these "short" distance scales: (1) what- 
ever is the technics used, recent works seem in agree- 
ment with "long" distance scales; (2) a difference in 
mass of 5 log M/Mq ~ —0.1 would mean that evolution- 
ary computations would b e completely wrong. We stress 
that iBeaulieu et al.1 l)200lh have the same point of view. 
Therefore we adopted a "long" distance value for the SMC 
distance modulus: /ismc = 18-9 ± 0.15 mag; this cho ice is 
supported by the recent result of H arries et al.1 l(2003h who 
found hsmg = 18.89±0.04 (statistical) ±0.15 (systematic) 
mag with a technics involving eclipsing binaries. 

Whereas the depth of the LMC seems to be negligible 
l|Van der Marel & Cionil[200l[k the depth of the SMC has 
been evaluated to range between ~ 0.2 and ~ 0.4 mag 
l|Crowl et al.ll200ll) . Then -for extreme cases- a given ob- 
ject inside the SMC could have an actual distance modulus 
+0.2 mag larger or lower than /zsmc = 18-9 mag, which 
has to be regarded as an average value. In order to esti- 
mate either the influence of an error on hsmc or an effect 
of SMC depth, we have made a test with fisMG = 19.1 
mag we got S\ogM/M Q ~ +0.1 and 5logL/L Q ~ +0.08. 

Another source of uncertainties is the reddening; if we 
assume an error of ±0.03 mag on E(B — V) = 0.08 mag, in 
turn we get a small uncertainty on masses and luminosity: 
<51og M/M e ~ ±0.01 and 5\ogL/L Q ~ ±0.05. 

Beside this, uncertainties connected to standardiza- 
tion of OGLE photometry are clearly negligible; with 
5V = ±0.02 mag we obtained 5 log M/M G ~ ±0.01 and 
SlogL/L G ~ +0.01. 

Moreover IBeaulieu et alJ (2001) made some additional 
tests: introducing turbulent convection, computing non- 
linear models or changing the meshes size within models 
does not yield to periods significantly different from those 
computed with LNA code. Therefore the uncertainty on 
distance (error on /zsmc or effect of SMC depth) remains 
the most important one. 



2.3.3. Comparison with ML-relations from evolutionary 
tracks 

From Fig. we remark a large discrepancy between the 
ML-relations derived from OGLE observations and from 
evolutionary calculations. For each evolutionary track, lu- 
minosity has been read at the "tip" of the blue loop, 
locus where the model spent at lot of t i me. T he dis- 
crepancy is also found using iBono et aD l)2000h mass- 
luminosity relation. The disagreement is getting worst 



when log M/Mq decreases. We must however emphasize 
that for log M/Mq ~ 0.4 (i.e. M/Mq ~ 2.5), the evo- 
lutionary track does not cross the Cepheid Instability 
Strip and a comparison between ML evo1 and ML puls for 
log M/Mq ~ 0.4 has no real meaning. Even an extreme 
value of /j,smc -he. 19.1 mag- can not lead to a perfect 
agreement between all ML-relations. 

3. Uncertainties in Standard Evolutionary Models 

In this section, we review the factors affecting the blue 
loop extension. Before presenting any models, we briefly 
recall a method allowing some predictio ns about the blue 
loop e xtension. We follow the work of lLauterborn et alJ 
ill 97ll) h who have defined an "effective core potential" : 



eff 



M, 
R 



£ (a Am AX) 



(6) 



where M c and R c are respectively the mass and the radius 
of the H c core and a a constant. Am is the width of the 
zone located between the H e core and the beginning of 
the outer chemically homogeneous region. AX represents 
the total hydrogen mass fract ion variation within Am . 
Numerical experiments done bv lLauterborn et alJ 1 197lil b 
have shown that a model undergoes a blue loop if this 



potential is lower than a critical value $ 



(crit) 
off 



We have to 



keep in mind this simple result: the lower $ e ff, the bluer 
the blue loop tip. 

In the next sections, we focus on a 3Mq track be- 
cause the most severe discrepancy in the CM-diagram is 
observed around this mass. 

3.1. Overshooting 

If we reduce the overshooting amount from a OVO r = 0.2 H p 
( "standard value" ) to a ovor = 0.0 H p , the H e core mass M c 
decreases as a consequence of the less extended H-core on 
main sequence. As a consequence, loops more extended to- 
ward the blue are expected. This is confirmed in Fig.^fb) 
where it clearly appears that even without any overshoot- 
ing (ao = 0.0 Hp), a 3 solar masses loop still remains too 
short to account for the observational data. 

3.2. Mixing Length Parameter 



The Mixing Length parameter «mlt = ^mlt /H p has been 
so far set equal to 1.6 in our stan dard models. Th i s valu e 
is derived from solar calibration l|Lebreton et 

ID £999)) 

and it is probably not universal: it may depend on metal- 
licity, mass, etc ... A priori, cumlt acts only on the con- 
vective flux (q?mlt is involved in MLT temperature gra- 
dient calculation) and does not change the position of 
Schwarzschild limit, hence AX in Eq. should remain 
unchanged, more blueward loops are not expected. The 
tracks computed with the extreme values «mlt = 1.0 and 
ctMLT = 2.0 are plotted in Fig. E^a). Both tracks have 
been calculated with a OVC r = 0.0 which is the most favor- 
able situation as explained in Sect. 13. II As one can notice 
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omlt has a negligible influence on the blue tip position. 
The effective temperature of the bluest point of the loop 
remains approximatively equal to ~ 3.76 (in Log) which 
is not enough to reach the warmer edge of the observa- 
tional IS, it still lacks ~ 850 K. We point out that a value 
of ttMLT = 1.0 is very unlike because it leads to a gi- 
ant branch around (V — I) ~ 1.5 where there are no stars 
within CM-diagram. 

3.3. Convective penetration 

Similarly, although not identically as the overshooting 
process, turbulent eddies must penetrate to some extent 
downward the convective envelope into stable radiative 
regions. However, we do not know how far they penetrate. 

Here we have carried out a calculation setting the ex- 
tension of convective penetrat ion at a cp = 0.7 H p fol- 
lowing the lAlongi et alJ l|199l|) prescription. They found 
that this value is needed to reproduce the properties of 
red giants branch luminosity function. This amount (0.7 
H n ) must be unde rstood as an order of magnitude as 
Alongi et al.l l)l99l|) 's calculations were performed before 
1992 when OPAL group published his new opacity tables. 

During the giant branch (hereafter GB) episode, the 
convective penetration produces a deeper penetration of 
the external convective zone. In this way, AX in Eq. [S] 
decreases and yields a lower $ c ft and bluer loop tip. 
Evolutionary tracks are displayed in Fig. E{c) for a 3 
Mq without overshooting (i.e. both with a ovcr = 0.0 
Hp) and show -as expected- that convective penetration 
slightly extends the loop but not enough to cross the en- 
tire IS. The extension difference reaches only a few ~ 130 
K, remaining too cold by ~ 720 K. 

3.4. Rotation 

iMaeder & Mevnetl (|200lh present evolutionary tracks in- 
cluding effect of stellar rotation at low metallicity Zq = 
0.004 in the mass range 9.0 to 60.0 Mq. The smaller 
mass value remains in t he Cepheid domain. In Fig. 6 of 
IMaeder fc Mevnetl <)200l[) the reader can remark that the 
blue loop extension is substantially reduced to logT e ff ~ 

4.12 to logics - 3.95 (AlogT cff 0.17) whereas we 

need AlogT c ff ~ +0.05. This blue loop reduction is due 
to the core extra-mixing added by convection equivalent 
to an overshooting addition. 

3.5. Influence of the CNO-cycle energy generation rate 

Although the CNO nuclear reactions cycle is rather well 
known, we have performed evolutionary calculations with 
an energy generation rate e cno artificially reduced by a 
factor two, from an arbitrary chosen post-main sequence 
stage (indicated by an arrow in Fig. |3|d). This magni- 
tude of uncertainty (a factor of two) is extremely large 
because the consequences during the Main Sequence phase 
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Fig. 4. Grid of evolutionary tracks for Zq = 0.001 with 
mass ranges between 2.5 M and 10.0 Mq. All the tracks 
cross the observational Instability Strip even the low mass 
tracks. Fundamental pulsators are represented by open cir- 
cles and first overtone ones by open triangles. 



would bring unavoidable disagreements between observa- 
tions and models. 

One can again predict what can be expected from such 
a numerical experiment. A lower e cno leads to a lower M c 
in Eq. El thus to a lower $ c ff and consequently to a bluer 
loop. Fig. |SJd) confirms this argument. Again, the loop is 
not extended enough and even an unrealistic uncertainty 
of a factor of two on the global energy generation rate e cno 
cannot explain the disagreement between observations and 
theory. 

Enhancing by a factor two the 3a reaction rate, is also 
found to have a negligible influence on the blue loop ex- 
tension. 
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3.6. Effects of helium 

The initial helium content adopted, Yq is expected to have 
only a minor influence on the blue loop extension, indeed: 

— the central helium content Y c during the blue loop 
episode does not depend on the helium content of the 
initial homogeneous model Yq. Hence Yq does not in- 
fluence £He the He-burning energy production rate be- 
cause within the inner regions Y c = 1 — Z c (with Z c the 
central heavy elements mass fraction) whatever the Yq 
value is. 

— during the blue loop the H-burning shell moves 
through the "A-profile" where Y varies between ~ 1 
(boundary of He core) and Y = Yq (chemically ho- 
mogeneous region mixed during the dredge-up episode 
when the model is closed to the Red Giants Branch). 
These Intermediate Y values are independent of Yq 
(obviously excepted values closed to Yo itself) 

Therefore the influence of Yo on the blue loop extension 
is expected to be very small. As a verification, models 
have been calculated with Yq = 0.23 and Yq = 0.28 which 
represent two extreme values: Yq = 0.23 is a rather low 
value for primordial helium and Yq = 0.28 which implies 
AYq/AZq ~ 9 while "r e asona ble" values are around 2, for 
a review see iLuridianal (|2002h . The tracks with Yq = 0.23 
and Yq = 0.28 in (Fig. Efe)) show that the initial helium 
content has no influence on the blue loop extension: the 
effective temperature of the tip remains equal to 3.76, i.e. 
~ 850 K colder than blue edge of IS, even in the favorable 
scheme of zero overshooting. 

3. 7. Effects of metallicity 

The high sensivity of a blue loop extension to metallicity 
is well known. The physical origin of this phenomenon is 
in the H-burning shell where material is processed through 
CNO cycle. For a fixed heavy elements mixture (here 
GN93) the lower Zq, the lower Ac, An, X<j (respectively 
C, N and O mass fractions) are. These three elements play 
the same role of catalysts in chemical reactions, there- 
fore a C, N, O deficiency leads to a lower energy gener- 
ated. Then, M c in Eq. © remains lower for a longer time 
and one obtains more extended blue loops. On one hand, 
the H-burning shell drives the star structure on the Giant 
Branch, on the other hand, during the blue loop episode, 
the He-burning core pulls the model towards the blue edge, 
where the He main sequence is located. The lower ccnO; 
the stronger He central burning effect. 

As a confirmation of this high metallicity sensivity, we 
have computed an evolution at 3 M taking a very low 
value: i.e. Zq = 0.001 which corresponds to [Fe/H] ~ —1.3. 
We compare the resulting extensions in Fig.^f). The blue 
loop crosses the entire IS, the tip reaching a position bluer 
than the blue edge of IS. 

Therefore the only way we have found to extend blue 
loops towards the high temperature edge of the HR dia- 
gram is to decrease the metallicity. In the next section we 



compare observational constraints and models built with 

Zq = 0.001. 



4. Models with Z = 0.001 




0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 

Log M/M Log M/M 



Fig. 5. Mass-Luminosity relations derived from OGLE 
observations assuming Zq = 0.001 corresponding to 
[Fe/H] ~ —1.3. Solid line: mass-luminosity relation from 
our evolutionary code with Zq = 0.001 and an overshoot- 
ing amount a OVC r = 0.2H p ; dot-dashed line: ML cvo1 for 
Zq = 0.004 and a ovcr = 0.2H p (same as Fig. [21 for com- 
parison). Panel (a): ML puls has been computed assuming a 
distance modulus ^smc = 18.9 mag; panel (b): same thing 
assuming (J-smc = 19.0. The diamond symbol shows the 
position of a 6 M Q , Zq — 0.001 model with a OVCI — 0.6 H p . 



4.1. Blue loops at Z = 0.001 

We have calculated a grid of evolutionary tracks at very 
low metallicity, i.e. Zq = 0.001. The results are dis- 
played in Fig. 01 where one can remark that the whole ob- 
served Instability Strip is crossed by the theoretical tracks, 
even the fainter part, i.e. the lower region of the color- 
magnitude diagram. These results suggest that a great 
part of SMC Cepheids could be metal deficient compared 
to the mean metallicity of the Small Cloud. 

One interesting point is that the shape of the 
Instability Strip at high magnitude is well reproduced by 
the decrease of the blue loop extension when going from 
higher to lower masses. 
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Fig. 6. Circles: Cepheids detected by OGLE team towards 
IC1613, evolutionary tracks (solid and dotted lines) has 
been calculated with Zq = 0.001 and an overshooting 
amount awr = 0.2 H p . 



4.2. Mass-Luminosity relation at Z = 0.001 

In Fig. 03 we have displayed the Mass-Luminosity rela- 
tions derived from OGLE observations assuming a metal- 
licity of Zq = 0.001. Fig. [SJd shows a better agreement 
between ML evo1 and MLP uls , if we assume msmc = 19.0 
mag -consistently with recent determinations- the agree- 
ment for low masses is excellent. Unfortunately it remains 
a discrepancy for higher masses -i.e. for log M/Mq ~ 0.7, 
this point will be discussed within the next section. 



5. Conclusion and discussion 

Sect. 0] shows that looking for an agreement between 
models and SMC observations for both blue loop ex- 
tensions and M-L relations for the SMC gives strong 
hints that high magnitude (i.e. low mass) SMC Cepheids 
could be metal deficient compared to the mean metallic- 
ity of the SMC; this fact could be explained by a "selec- 
tion effect" : only stars with low enough metallicity could 
have an evolutionary track crossing Cepheids Instability 
Strip. Unfortunately direct spectroscopic determinations 
of [Fe/H] for the SMC Cepheids around I-mag ~ 17 are 
not yet available. We discuss now a few issues in favor 
of the above proposition. First we consider another metal 
poor galaxy and then discuss the information which can be 
drawn from Cepheids in SMC clusters and beat Cepheids 



which supports the existence of SMC Cepheids as metal 
poor as [Fe/H] ~ — 1.3 dex. 

5.1. Comparison between our model and a very low 
metallicity Cepheid population 



Udalski et al provided a sample of Cepheids be- 

longing to the galaxy IC1613 w ith a metallicity [Fe/H] ~ 
-1.3 ± 0.2 fsee lLee et al"1ll993j) . This data set offers the 
opportunity to check whether our evolutionary models 
are valid for metallicity as low as the value suspected 
for SMC Cepheids located below the A line in a CM- 
diagram. Fig. displays OGLE objects and our evolu- 
tionary tracks for Zq = 0.001. The distance modulus for 
IC1613 is tak en to be 24.2 ± 0.1 mag and absorption is 
A| = 0.05 (seelUdalski et al.N2001l) : reddening is given by 
ISchlegel et al.l l)l998|) . Fig. shows that the evolutionary 
tracks cross the whole instability strip as defined by the 
observed Cepheids. For this galaxy, no problem of blue 
loop exists with our models indicating that the main fea- 
tures of our models likely capture the essential evolution- 
ary properties at low metallicity. One can notice that even 
a cut-off of OGLE detector around I-mag ~ 22 can not re- 
ally change our conclusion because the blue loop for 2.5 
Mq is extended enough to model Cepheids as faint as I- 
mag ~ 22.2. 

5.2. Evidences for metal poor stars within the SMC 



iLuck et all l|l998f) have determined the chemical composi- 
tion of Cepheids in SMC by means of high resolution spec- 
troscopy. The [Fe/H] values found by the authors range 
between —0.84 and —0.65 corresponding to a mass frac- 
tion Zo ranging between ~ 0.0030 and ~ 0.0045, these 
values bracket the commonly assumed metallicity mean 
value for SMC, i.e. Zq — 0.004. Howev er two rema r ks are 
in order here: (1) the sample studied bv lLuck et all lll998h 
has a quite poor statistic (6 objects) while OGLE sam- 
ple contains about 2000 objects; (2) m ore importantly, as 
shown by Fig.^b) the stars studied bv lLuck et al.l 1 1998) 
are objects much brighter than those around I-magn^ 17 
which are not reached by our blue loops. The reason for 
the choice of bright objects is that spectroscopic determi- 
nations are easier for brighter objects. Therefore any biais 
-concerning high magnitude SMC Cepheids metallicity- 
cannot be excluded. 

On the other hand, one may think that it is possi- 
ble to infer some indirect information about the SMC 
Cepheid metallicity. OGLE team has indeed discovered 
many stellar clu sters in the SMC, (see IPietrz vnski et all 
119981) . Moreover IPietrzvnski fc Udalski ljl999d) have de- 
tected 132 Cepheids belonging to these clusters. One of 
the main properties of stars belonging to a given cluster is 
to present the same chemical composition. Therefore any 
indications about metallicity of these SMC clusters gives 
an information about the metallicity of Cepheids belong- 
ing to clusters. The literature is quite poor about metal- 
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licity deter minations for SMC clusters. Table ^ mainly 
taken from ICrowl et al.l l)200l|) gives metallicity estima- 
tions for SMC clusters. Thanks to a cross identification 
we have found 2 clusters belo nging to the catalogue of 
IPietrzvnski fc Udalskil l|l999cfl and having a metallicity 
determination in the literature: NGC 330 (SMC0107 in 
OGLE catalogue) and NGC 416 (SMC0158). However these 
clusters are known to be young or intermediate-age sys- 
tems; one has to check whether the ages of Cepheids are 
compatible with age of the cluster hosting them. 

In the case of NGC 330 IChiosi et all l|l995[) 
have found a maximum age of 48 Myr; beside this 
IPietrzvnski fc Udalskil 1 1999b) have deriv ed from their 
study 3 1.6 + «k Myr. From HST observations lMighell et"all 
1 1998b() estimate an absolute age for NGC 416 of 6.6 ± 
0.5 Gyr assuming that L i ndsay 1 cluster is 9 Gyr 
old; IPietrzvnski fc Udalskil l|l999b|) confirm that NGC 
416 is older than 1 Gyr. We underline that neither 
NGC 330 nor NGC 416 are mentioned in the erratum 
IPietrzvnski fc Udalskil l|l 999ah . 

For the Cepheid (SMC_SC7 206038 in OGLE cata- 
logue) suspected to belong to NGC 330 if we assume a 
metallicity [Fc/H] = —0.7 (Z = 0.004) and a mass about 
3.5 M Q consistently with its position within CMD; we 
found an age of ~ 230 Myr; assuming [Fe/H] = — 1.3 
[Z = 0.001) and 3.0 M Q we obtain - 300 Myr. Thus 
the Cepheid is likely a field star and does not belong to 
NGC 330. 

Stars suspected to be NGC 416 objects have an I- 
magnitude between 15 and 16 mag; this corresponds to a 
mass around 4 M for [Fe/H] = —1.3 (metallicity consis- 
tent with cluster one). From our evolutionary calculations 
we get -for such mass and metallicity- an age of ~ 150 
Myr. 

Unfortunately this value is not compatible with the es- 
timated age for NGC416; indeed the age spread (probably 
around ~ 0.5 Gyr) does not allow such young objects to 
belong to the cluster. Nevertheless one can notice that the 
majority of SMC clusters in Tab. are metal deficient, 
thus SMC stars with metallicity lower than Zo ~ 0.004 
([Fe/H] ~ —0.7) does exist and the hypothesis of metal 
deficient SMC Cepheids appears to be reasonable. 

5.3. Information brought by SMC Beat Cepheids 

OGLE Tea m has discovered a sa mple of 93 beat Cepheids 
in the SMC. IUdalski et all l|l999h found that 23 pulsate si- 
multaneously on the fundamental mode (hereafter F) and 
the first overtone (hereafter 10T), the remaining objects 
have been found to pulsate simultaneously on the first and 
second overtones (hereafter 20T). 

In order to derive their ML puls iBeaulieu etaD l|200lj) 
choose three quantities among the four observational ones: 
T c ff, L and the periods P k and P k+1 (P for F/IOT 
and Pi for 10T/20T); they calculate the theoretical 
value of Pk+i which is noted Pk+i(calc) (the observed 
one being Pfe + i(obs)). They next define the parameter 



Table 1. Determinations of SMC clusters metallicity 



Cluster 


[Fe/H] 

L / J 


Cepheid(s)? 




a oo _i_ nil , ■. 
— U.OZ ± U.ll (a) 


no 


"\T /"I /"I A 1 1 


— U.OO ± VJ.Ul (b) 


no 


IN^jV-' ±DZ 


A OA _|_ n 1 ^ / \ 
— U.y4 it U.1D (c) 


no 


Lindsay 113 


1 oa xn 11 / n\ 
— 1.Z4 ± U.ll (d) 


no 


Kron 3 


1 1 d A AO / i \ 

— l.lO ± U.Uy (d) 


no 




— l.DU ± U.14 (d) 


no 


TV/" 1 A1 R 


1 A A -4- fl 10 / j\ 
— 1.44 ± U.1Z (d) 


no 


NGC 361 


-1.45 ±0.11 (d) 


no 


Lindsay 1 


-1.35 ±0.08 (d) 


no 


NGC 121 


-1.71 ±0.10 (d) 


no 


Kron 28 


-1.20 ± 0.13 (e) 


no 


Lindsay 38 


-1.65 ±0.12 (e) 


no 


Kron 44 


-1.10 ±0.11 (e) 


no 



(a'Hill 11999); (b Alvcs & Saraiodini (1999); (V i 
fd lMighell et 3,1,1115553) : Ce lPiatti et alj ISSCT) 
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Fig. 7. Numerical experiments for a 6 M© model. For both 
metallicities Zq = 0.004 (panel a) and Zo = 0.001 (panel 
b); we have computed a grid of evolutionary tracks varying 
the overshooting amount from 0.2H p up to 0.6H p . For high 
metallicity blue loop disappears suddenly for overshooting 
larger than 0.4H p while the blue loop extension decreases 
monotonously in the low metallicity case. 



e = Pfe+i(calc)/P/c+i(obs)) allowing a comparison be- 
tween theory and observations. They explore the influ- 
ence of different important parameters, particularly dis- 
tance modulus and reddening and finally conclude that a 
solution (i.e. e ~ 1) is found simultaneously for F/IOT 
pulsators and 10T/20T only if the metallicity is settled 

to Z Q = 0.001 (i.e. [Fe/H] 1.3). 

In Fig. [S| we have plotted the observed beat Cepheids 
together with evolutionary tracks and the straight line A 
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defined in Sect. [21 This plot clearly shows that the beat 
Cepheids are located in the region where we suspect that 
objects are metal deficient (i.e. with a metallicity around 
Zq ~ 0.001). All 10T/20T pulsators are below A (ex- 
cepted one object) while F/IOT pulsators are scattered 
slightly above and below. 

Thus these pulsation/evolution models of SMC beat 
Cepheids argue in favor of a relation between a metal de- 
ficiency (with respect to the mean value of the SMC) and 
the existence of SMC Cepheids at low magnitude. 

5.4. The case of high mass Cepheid 

Although it is slightly out of scope of this paper where 
we focus on faint SMC Cepheids, we will debate in this 
paragraph the case of brighter objects: i.e. log M/M Q ~ 
0.7- - 0.8 corresponding to 3.6 < logi/L < 4.1. 
Whatever is the assumed metallicity: Zq = 0.004 or 
Zq = 0.001, few objects with an evaluated mass around 
log M/Mq ~ 0.77 (M ~ 6 Mo) have a MLP uls in discrep- 
ancy with the related ML cvo1 . 

We can make some hypothesis: as shown in Sect. 12.3721 
with a larger distance modulus -i.e. ^smc — 19.1 mag- we 
get S\ogM/M Q ~ +0.1 and 5 log M/M Q ~ +0.08 (com- 
pared with the situation with 18.9 mag), but this extreme 
value is not able to bring a full agreement between ML puls 
and ML cvo1 for brighter objects. On the other hand it is 
unlike that all these stars would be located deeper in SMC 
than others. 

For a given mass value, a way to enhance the luminos- 
ity is to c onsider a larger over shooting: amount. As sug- 
gested by ICordier et al.1 l|2002j) . it can not be excluded 
that average overshooting amount for intermediate mass 
stars increases when metallicity decreases; they have de- 
rived -assuming Zq = 0.004 for SMC main sequence stars- 



0.40 



" nag H r . W ith a LMC bump Cepheids 



studv lKeller fc Woodl l|2002l) infer an overshooting amount 
A c = 0.63 ± 0.03 Hp (~ 0.3 H p in our formalism). 

We made a test -involving all mass values- with a OVO r = 
0.4 Hp, as expected blue loop extensions are reduced 
(for low masses tracks the excursion of blue loop within 
IS is less deep) and luminosity is increased not enough 
(SlogL/ 'Lq ~ +0.1) to get an agreement between ML cvo1 
and ML puls for log M/M ~ 0.77. Then, to increase the 
overshooting amount over the whole range of mass is not 
the solution. 

Another possibility is that overshooting can depend 
on mass, increas i ng wi th mass; as suggested for instance 
by lYoung et al ] ll200llL Thus we have concentrated us 
on 6Mq models, varying overshooting amount between 
0.2 Hp (our "standard" value here) up to 0.6 H p ; this for 
both metallicities: Z = 0.001 and Z = 0.004. Results 
are shown in Fig. [7| it is clear that beyond 0.4 H p one 
no longer gets blue loop for Zq = 0.004. In contrast with 
Zq = 0.001 blue loop extension decreases, but in a way 
where Cepheids IS is crossed by tracks. logL/L increases 
reaching large enough value (for 0.6 H p , see diamond sym- 



bol in Fig.EJb)) in order to make it up M L evQl and ML puls 
for lo g M/Mq ~ 0.77. We stress that iKeller fc Woodl 
(2002) found their quite "high" overshooting amount us- 
ing a sample of bright LMC Cepheids, this support our 
proposal of an higher overshooting for brighter Cepheids. 
One more time our work favors a solution involving low 
metallicity. Towards high masses, another selection effect 
could occur if the overshooting increases with mass and 
reached ~ 0.6 H p for masses larger than ~ 6M at low 
metallicity. A detailed study is needed on this topic and 
is out of the main goal of this paper. 



5.5. Summary 

In this paper we have explored two main problems related 
to the SMC Cepheid population: (1) the blue loop exten- 
sion for high magnitude stars, (2) the Mass-Luminosity re- 
lation. We have first shown that blue loop extension is ex- 
tended enough only if the metallicity is substantially lower 
than the commonly used value for SMC objects models. 

Evolutionary tracks computed with Zq = 0.001 cor- 
rectly reproduce the Instability Strip shape for low masses 
and Mass-Luminosity relation derived from these tracks is 
in rather good agreement with ML p uls deduced from obser- 
vations using a technique similar to lBeaulieu et alJ l|200ll) 
one. The remaining discrepancy for the small population 
of brighter objects could be explained by a joined effect 
of low metallicity and rather enhanced core mixing pro- 
cess. Further re searches are needed on this subject. We 
emphasize that IPietrukowicd l)2002|) -who estimates pe- 
riod change rates of SMC OGLE Cepheids- found also for 
brighter objects a rather bad agreement between models 
and observations. Our point is that all Cepheids below 
A are likely metal poor and stars above A belong prob- 
ably to a "mixed" population. Finally, the present work 
strongly suggests the existence of an evolutionary selection 
effect for fainter Cepheids belonging to SMC. High reso- 
lution spectroscopic chemical composition determinations 
for SMC Cepheids through the entire IS and particularly 
around magnitude 17 are requested to bring definitive ar- 
guments in favor or against the present suggestion. This 
could be possible for a sample of few stars with UVES 
VLT spectrograph. Results would put a new light on the 
cosmologically important metallicity dependence of the 
Cepheids Period-Magnitude relation. 
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Fig. 3. Influence of free parameters on the blue loop extension for a 3 MQinodel. The parameters are: (a) the mixing 
length, (b) the overshooting, (c) the convective penetration (overshooting below external convective zone), (d) CNO 
nuclear cycle energy production rate, (e) initial helium content and (f) initial metallicity. For each plot, a vertical 
segment (defined in Fig. ^p) shows the temperature which the loop must reach to cross over the entire observational 
IS. 



